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Correlations may affect propagation processes on complex networks. To analyze 
their effect, it is useful to build ensembles of networks constrained to have a given 
value of a structural measure, such as the degree-degree correlation r, being random 
in other aspects and preserving the degree sequence. This can be done through Monte 
Carlo optimization procedures. Meanwhile, when tuning r, other network proper- 
ties may concomitantly change. Then, in this work we analyze, for the r-ensembles, 
the impact of r on properties such as transitivity, branching and characteristic path 
lengths, that are relevant when investigating spreading phenomena on these net- 
works. The present analysis is performed for networks with degree distributions of 
two main types: either localized around a typical degree (with exponentially bounded 
asymptotic decay) or broadly distributed (with power-law decay). Size effects are 
also investigated. 
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I. 



INTRODUCTION 



Complex networks are realistic substrates for simulating many social and natural phe- 
nomena. To address the influence of network topology, primarily, different classes of degree 
distributions P(k) can be considered. Meanwhile, for a given distribution of degrees, cor- 
relations may give rise to important network structure effects on the studied process [THE]. 
These structural effects may have important consequences, for instance, correlations may 
shift the epidemic threshold [6]. Although correlation effects may be absent in some cases 
[7J, in other ones, they can not be neglected. 

Despite there are efficient algorithms to generate networks with fixed degree-degree corre- 
lations [8], real joint probabilities of two or more degrees measured in networks of moderate 
size may be noisy and hard to be modeled. Then, operationally, average nearest- neighbors 
degree distributions |9j or single quantity measures are used. Although other variants have 
been defined in the literature [TO], HTj . as quantifier of the tendency of adjacent vertices to 
have similar or dissimilar degrees, we will consider the standard measure of (linear) degree- 
degree correlations, namely, the assortativity (Pearson) coefficient [12] 



where (• • • ) e denotes average over edges and k and k' are the degrees of vertices at each 
end of an edge. Despite this coefficient is known to present some drawbacks [13j, it is a 
standard and commonly used quantity, hence being worth to be analyzed. Moreover, it has 
the advantage of being a single value measure, that is easier to be controlled than other 
multi-valued quantities. 

To analyze the influence of correlations, as well as of any other structural feature, it is 
useful to build ensembles of networks holding that property, while keeping fixed the sequence 
of degrees. As it will be described in Sec. [TTJ this kind of ensembles can be achieved by means 
of a suitable rewiring, performed through a standard simulated annealing Monte Carlo (MC) 
procedure to minimize a given energy-like quantity (maximum entropy ensemble approach), 
function of the graph property to be controlled (r in our case) [T1HT6] . Once tuned r, it 
is important to characterize how other network properties are altered as by product. Some 
interdependencies among certain network properties have already been numerically shown in 
the literature, for real as well as for artificial graphs [15]. Analytical relations have also been 
derived [TTHT9] . Because of its crucial role in spreading phenomena [20], we will focus here 



r = 



(kk% ~ (k) 2 e 

(k 2 )e ~ (k)l ' 



(1) 
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on the effect of r over typical distance measures as well as on the branching and transitivity 
of links. 

As a measure of the average separation between nodes, we consider the average path 
length [21]. In the subsequent calculations we use the expression, 

T = £r =1 (£i)aw-i) (9) 

where n is the number of (disconnected) clusters and Ni is the number of nodes in cluster i. 
Moreover, being dkj the distance (number of edges along the shortest path) between nodes 
k and j (taking d^j = if the nodes do not belong to the same cluster), then 

Alternatively, in order to avoid the issue of the divergence of the distance between dis- 
connected nodes, we consider the inverse, 1/E, of the so-called efficiency [22] 

where N is the number of nodes. It represents a harmonic mean instead of the arithmetic 
one. We also compute the diameter D = max{dy}. 

The transitivity of links can be measured by the clustering coefficient [23], |2l] 

C=^^ , (5) 

J2i=i ki(ki 1) 

where ha is the number of triangles and ki is the degree of node i. We also considered the 
mean value, C, of the local clustering coefficient Cj, defined as Cj = 2ei/{ki{ki — 1)), where 
Ci is the number of connections between the neighbors of vertex i [21 J. We took Cj = when 
ki — or 1. Other measures that arise in the decomposition of r [17J will also be considered. 

Besides detecting interdependencies among structural properties, it is also important to 
know how these properties depend on the system size N. We will analyze these issues for two 
main classes of degree distributions (Poisson and power-law tailed). We will also investigate 
real networks degree sequences. 



II. NETWORKS AND ENSEMBLES 



For each class of networks, we will consider different values of the size, N, and the mean 
degree, (k), within realistic ranges. 



4 



As a paradigm of the class of networks with a peaked distribution of degrees, with all 
its moments finite, we consider the random network of Erdos and Renyi [2H]. Within this 
model, a network with iV nodes is assembled by selecting M different pairs of nodes at 
random and linking each pair. The resulting distribution of links is the Poisson distribution 
P(k) = (k) k /k\, where the mean degree is (k) = 2M/N. 

We also analyze networks of the power-law type, i.e., with P(k) ~ /c -7 , 7 > 2, corre- 
sponding to a wide distribution of degrees, with power-law tails. Then, moments of order 
n > 7 — 1 are divergent. We built power-law networks by means of the configuration model 
|26j . Following this procedure, one starts by choosing iV random numbers k, drawn from 
the degree distribution P{k). They represent the number of edges coming out from each 
node, where these edges have one end attached to the node and another still open. Second, 
two open ends are randomly chosen and connected such that, although multiple connections 
are allowed, self connections are not. This second stage is repeated until each node attains 
the connectivity attributed in the first step. If eventually an edge has an open end, then it 
is discarded. However, for large networks, the fraction of discarded edges is negligible. To 
draw the set of numbers k with probability P{k) = A/7c~ 7 , with k m i n < k < k max (hence the 
normalization factor is Af = 1/Ylk max ^~ 7 )' we used the inverse transform algorithm [27]. 
Notice that k max < N — 1 and k max » k min , then we determined k m i n to fit the selected 
value of (k) (within a tolerance of at most 1%), such that 

/U\ — ^ k min ^ I t ""max ^min ^ J_ L _i ,c\ 

Efc:; fc-7 7-2 kiai - k i m p n 7 - 2 

It is worth mentioning that the value k = N — 1 is not usually achieved, the natural cut-off 
being k c ~ [28] . 

In order to attain a desired value of r, we follow an standard rewiring approach. We 
want to build an ensemble of networks {G} with a given value of r (r-ensemble) but that are 
maximally random in other aspects, i.e., making the fewer number of assumptions as possible 
about the distribution P(G). Then, we use an exponential random graph model, such that 
the set of networks {G} has distribution P{G) oc e~ H ^ G \ where H(G) is a Hamiltonian or 
energy-like quantity [l4j. In order to get an r-ensemble, with r = ?\, we consider [15] 



H(G) = P\r-n\, (7) 
where /3 is a real parameter. The ensemble can be simulated by means of a MC proce- 
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dure: at each step, a rewiring attempt is accepted with probability min{l, e~^ H ( G ) _i? ( G )l}. 
Rewiring steps are performed by randomly selecting two edges that connect the vertices 
a, b and c, d, respectively, and substituting those two links by new ones connecting a, c 
and b, d [29] . Movements yielding double links are forbidden. Notice that this process 
preserves the connectivity of each node. We start the simulation by taking f3 = [during 
at most 100 MC steps (MCS), where each MCS corresponds to N attempts]. The effect 
of this stage is basically to destroy multiple edges. We did not notice any clear hysteresis 
effect like those observed when controlling, instead, the number of triangles with a differente 
Hamiltonian [30]. Subsequently, is increased (in increments A/3 = 1000), at each 50 MCS, 
until r stabilizes, typically attaining the prescribed value r*. Then, the quantities of inter- 
est are calculated and the whole process is repeated, starting with a new degree sequence. 
For power-law degree distributions, we observed that the process is non-ergodic, hence we 
computed sample mean and standard deviation over 100 realizations of the described pro- 
tocol. We checked that the choice of other expressions for H(G), vanishing at r*, did not 
significantly affect the results but just the convergence time. 



III. RESULTS 

Let us start by reporting the effects of r on the clustering coefficient C. 

For the Poisson case, we depict in Fig. [IJa) the behavior of C as a function of r for a fixed 
number of nodes (N = 8000) and different values of the mean connectivity (k). Very small 
values of C emerge. The transitivity C monotonically increases with r. This is consistent 
with the results of Ref. [15] (restricted to r > 0) for such kind of networks. We observe two 
regimes with a crossover at r ~ 0.5: a very slight increase with r below the crossover and a 
more pronounced one in the region above it. The existence of two regimes could be related 
to the assymetric character of r, which does not measure assortativity and disassortativity 
on the same grounds. Below the crossover, C linearly increases with (k) about one order 
of magnitude within the analyzed range. Meanwhile, above the crossover, C remains of the 
same order when the average connectivity increases, even for small (k) (also see the inset of 
Fig. [T|a) where C is plotted vs (k) for selected values of r). In Sec. [TV] we will discuss these 
issues in more detail. For the mean local clustering coefficient C, we obtained a qualitatively 
similar dependence on r than that observed for the clustering coefficient C . However, the 
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FIG. 1. Clustering coefficient C as a function of r for Poisson networks: (a) N = 8000 and 
different values of (k) indicated on the figure. The graph shows a monotonic increase of C with r. 
There are two regimes: a very slight increase of C with r below r ~ 0.5 and a more pronounced one 
above it. (b) (k) = 4 and different number of nodes N, also indicated on the figure. As the number 
of nodes increases, C decays with the asymptotic law C ~ 1/N, characteristic of uncorrelated 
random graphs. Standard errors are about 10%. Dotted lines are a guide to the eyes. The insets 
show C vs (k) (a) and N (b) for selected values of r (-0.6, 0.0, 0.6 and 0.8). 



increase of C with (k) is linear for any fixed r. For r = 0, C = C = (k)/N, as expected. 

In Fig. []Jb), size effects are exhibited for (k) = 4, representative of the other values 
considered. As the number of nodes increases, C decays as C ~ 1/N for all r (as depicted in 
the inset). Therefore, in the r-ensemble of Poisson networks, transitivity is only a finite-size 
effect and vanishes in the infinite network (thermodynamic) limit with the same asymptotic 
law C ~ 1/N that for an uncorrelated random graph [3 lj . 

For the power-law class, the range of allowed values of r is restricted. That is, values 
of r arbitrarily different from zero can not be attained in typical realizations of the MC 
protocol described in Sec. |nj In order to determine the typical maximal (minimal) values, 
r max ( r min), w e imposed 7* = 1 (-1) and detected the stationary values of r. The time 
evolution of r for r ± = 1 (-1) is illustrated in Fig. |2](a) for 7 = 3.5, = 5000 and (k) ~ 4. 
Notice the large deviations amongst the steady values of different realizations mainly for 
the upper bound. We verified that this picture does not change by implementing other 
definitions of H(G) in Eq. Q, e.g., (3\r - r*\ a , with a ^ 1. Average extreme values (over 
100 samples, after 3 x 10 4 MCs) are displayed in Fig. [21(b) , as a function of the system size, 
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for different values of 7. For fixed size, the lower 7, the narrower the allowed interval of 
r. In fact, in networks constrained to a given degree sequence, structural limitations (or 
correlations) arise: either multiple connections or dissortative two vertices correlations [32J. 
For instance, the exclusion of multiple connections hampers the natural tendency that hubs 
connect among them, hence diminishing the assort at ivity. This effect is more pronounced 
the smallest 7. For fixed 7, the interval shrinks with system size, for 7 < 3, due to the 
divergence of fluctuations in the large N limit [13] - In Ref. similar restriction was also 
observed for 2 < 7 < 3, although instead of the average connectivity, k min was kept constant 
(kmin = 6). In that case, it was reported that the upper and lower bounds both tend to 
zero, hence r — > in the infinite network limit. In fact, we observe that in that interval of 
7 (e.g., 7 = 2.5) both bounds are negative and as N increases the allowed interval collapses 
to a negative value that tends to zero. We noticed restriction in the correlation bounds for 
7 > 3 too. As N grows, the lower bound also increases towards zero or at least to a small 
finite value. Simultaneously, the upper bound seems more stable, however its asymptotic 
behavior is not neat yet, even having considered up to iV > 10 5 . Moreover, as N increases, 
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FIG. 2. The range of allowed values of r is restricted for the power-law class: (a) Time evolution of 
r, after setting r* = 1 (-1) to obtain r max (r m j„), for networks with power-law degree distribution 
(with 7 = 3.5 and iV = 5000). Shown are 8 individual samples (thin lines) and their respective 
averages (thick lines) . (b) Average extreme values [r max (open symbols) and r m in (filled symbols) , 
standard errors are at most 50%] vs system size iV for different values of 7 indicated on the 
figure. Dotted and dashed lines are guides to the eye for r max and r m j n , respectively. In all cases 
(k) = 4.00 ± 0.04. For a given size, the allowed interval of r is narrower for lower 7. 
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it takes longer to attain steady states. 

The allowed interval of r is quite restricted for scale-free networks. However, we still 
analyzed systematically cases with 7 > 3 (7 = 3.5, 4.0 and 4.5), yielding finite second 
moment. Even, in this range, the accessible interval of r is limited, then, we proceeded 
as follows. If the desired r* is not attained, within a tolerance of 10 -3 , in 2 x 10 4 MCS, 
that instance is discarded and we make a new trial. If we did not attain 100 successes in 
200 trials, the procedure is interrupted. Alternatively to the configuration model, we also 
started from networks generated by preferential attachment [33], yielding similar results. 
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FIG. 3. Clustering coefficient C as a function of r, as in Fig. [T] but for power-law networks with 
7 = 4.0: (a) N = 8000 and different values of (k). In the assortative region, C reaches larger values 
than in the Poisson class, (b) (k) = 4 and different number of nodes N. Note that for assortative 
networks a finite degree of clustering seems to persist for large networks. Standard errors reach 
50% for the lowest values of C. In the inset, missing values are due to the limitation in attaining 
the prescribed values of r. 



The outcomes for the power-law class with 7 = 4.0 are displayed in Fig. [3] Standard errors 
are larger in the power-law case, likely due to the variability in the tails of the distribution 
of links from sample to sample. Outcomes for the other two values of 7 studied (3.5 and 
4.5) display features similar to those of the case 7 = 4.0 used as illustrative example, despite 
the third moment becoming divergent at 7 = 4.0. Two regimes are also observed, with the 
crossover now closer to r = 0, but some qualitative differences appear in comparison to the 
Poisson case. C rapidly increases with r, attaining, for assortative networks, larger values 
than in the Poisson class. These large values are due to the inclusion of highly connected 
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nodes, absent in the Poisson networks, that for large r* tend to gain links among them, 
contributing strongly to r and also to C . 

With respect to finite size effects, below the crossover the small non-null C is again 
due only to the finite size of the network. However, in the assortative region (above the 
crossover), it seems that a finite degree of clustering persists for large networks (see inset of 
Fig. gb)), in contrast to the Poisson case and to the dissassortative region. In fact, notice 
that, when iV increases one order of magnitude, C decreases also one order of magnitude 
in the dissortative region, while C remains of the same order in the assortative interval. 
Even if C vanished in the infinite size limit, since the decay is very slow, then an effective 
clustering would remain in moderate, or even large, size networks. We will discuss the 



interplay between C and r further in Sec. IV For the mean local clustering coefficient C, 
a qualitatively similar dependence on r is observed, but with smaller values. Moreover, C 
increases linearly with (k), in the analyzed range, for any fixed r, not only for dissortative 
networks, and C decays with N for any r. Forr = 0, C — C — ({k 2 } — (k)) 2 / (N (k) 3 ) [T9] . 
as expected. 

Let us analyze now the influence of r on network characteristic lenghts. The depen- 
dency of the measures 1/E, L and D on r is depicted in Fig. |4j for Poisson and power-law 
distributed networks, with N = 8000 and < k >~ 4. 1/E and L have close values, system- 
atically shifted. In first approximation both types of network yield similar values of 1/E 
(hence also L), given N and (k). However, the diameter D is more dependent on the type 
of network. It is larger and is more strongly affected by r in the homogeneous Poisson case. 

Plots of 1/E vs r for different values of N and (k) are shown in Figs. [5] and [6] for Poisson 
and power-law networks respectively. In both cases, the networks display the small-world 
property [2T] (even smaller in the power-law case) with a slow (logarithmic) increase with 
N and a smooth decrease with (k) (see insets of Figs. [5] and [6]). 

In the Poisson case, the mean path does not depend on r significantly when (k) is not too 
small (> 6), as indicated by the relatively flat plots in Fig. |5](a). Only for small (k) there 
are important effects for assortative correlations. For instance, for (k) ~ 4 (Fig. |5](b)), 1/E 
increases in about two units from r ~ to r ~ 1, for all the analyzed sizes. This effect is 
still larger for (k) = 2 where L increases by a factor about two in the same interval of r, as 
shown in Fig. ga) for N = 8000. 

In order to further interpret these results, we investigated the cluster structure of the 
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FIG. 4. Distance measures 1/E, L and diameter D, as a function of r for Poisson (open symbols) 
and 7 = 4.0 power-law (filled symbols) networks. In all cases N = 8000 and (k) = 4.00 ±0.04. The 
inset is a zoom of the main plot. At first sight, the two types of network display similar values of 
1/E and L, for a given N and (k). The diameter is more sensitive to the type of network and is 
more affected by r in the Poisson case. 
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FIG. 5. Mean distance 1/E as a function of r for Poisson networks: (a) N = 8000 and different 
values of (k) indicated on the figure. The networks exhibit the small-world property, (b) (k) = 
4.00 ± 0.04 and different number of nodes indicated on the figure. The effects of r on the mean 
path are significant only for small (k) because of the fragmentation of the network. 



resulting rewired networks. We measured the size of the largest cluster (let us call it N±), 
the number n of clusters, and the average size S of the clusters different from the largest 
one. The plots are presented in Fig. [7] for (k) = 2 and 4. For (k) = 4, the relative size 
of the largest cluster (circles) is about 0.98 for most of the range of r, notice however that 
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FIG. 6. Mean distance 1/E as a function of r as in Fig. [5] but for power-law network with 7 = 4.0. 
In this case, paths are shorter than in the Poisson nets. 



it slightly decays towards \r\ = 1 (which is more evident for (k) = 2). As (k) increases, 
the number of fragments rapidly decays and the average size S (triangles) tends to unit, 
meaning that only single nodes are disconnected (recall also that P(0) = e - ^). Therefore 
the increase of the mean distance towards |r| = 1, observed for low (k) in Fig. |5j may simply 
reflect the fragmentation of the network. Clearly, for high values of the assortativity, the 
network tends to fragment into groups of vertices that have the same degree. Moreover, 
for values of (k) approximately larger than 1, the percolation analysis performed by Noh 
on Poisson networks [16] shows that the size of the largest cluster is smaller for assortative 
networks than for dissortative and neutral ones. Meanwhile, as (k) increases, the fraction of 
vertices that do not belong to the largest cluster becomes negligible, although more slowly 
the more assortative the network. Therefore, in such large (k) limit, the mean distance 
remains invariant under changes of r. Hence, transport processes modeled in these networks 
may suffer important impact when r is large and (k) small. The longer the typical separation 
between nodes, the slower the propagation. 

In power-law networks, for fixed N and (k) > 2, there is an interval of r where paths 
are shorter than in the Poisson nets (Fig. [4]), and still shorter as 7 decreases (not shown). 
Moreover 1/E becomes more sensitive to the coefficient r (smile shape), in the region where 
plots are flat for Poisson networks. Notice also that minimal mean paths occur for slightly 
assortative correlations (r > 0), slowly increasing with iV (Fig. [6^b)). The analysis of 
clusters for 7 = 4.0, shows that for (k) > 4 there is a single cluster of size N, for all r. Only 
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FIG. 7. Clusters analysis. Plots of N\/N (where Aq is the size of the largest cluster) [circles], 
average size of finite size clusters, S [squares], and number of clusters, n [triangles], as a function 
of r for Poisson networks with (k) = 2 (open symbols) and 4 (filled symbols). The fragmentation 
of the network observed for high values of the assortativity reflects the grouping of vertices with 
the same degree. The figure shows outcomes for N = 8000. Outcomes for sizes N = 4000, 8000 
and 16000 all collapse into single curves (not shown). 

for (k) = 2.0 we observed fragmentation with Aq/iV ~ 0.7 — 0.8, n/N ~ 0.06, S ~ 4 for all 
N > 2000 (plots not shown). 

For the mean path L, we observed qualitatively similar outcomes although shifted to 
slightly higher values, as illustrated in Fig. |4j 

We also applied the rewiring procedure described in Sec. [TT] to real degree sequences. 
Networks were symmetrized and edge weights were ignored. In Fig. [8] we depict the behavior 
of C and 1/E vs r, for the PGP (encrypted communication network using Pretty Good 
Privacy encryption algorithm) largest component [34], P2P (Gnutella peer-to-peer network) 
[35] . PGR (electrical power grid of the western United States), example of small- world 
network [36] and APC (astrophysics collaboration network) [31]. First notice that in all cases 
the clustering C is much larger in real networks than in the randomized (r-ensemble) ones, as 
already observed for other examples in Ref . [Po] . The mean distance is also typically larger 
in the real networks. An exception is P2P network, characterized by a value of 1/E typical 
of the r-ensemble. Rewired real networks display some of the typical behaviors observed for 
the artificial cases. Let us make some remarks arising from comparisons, (i) PGR (power 
grid) displays plots of C vs r and 1/E vs r that are in good accord with those observed 
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FIG. 8. C and 1/E vs r, for real networks. Original values, before rewiring, are also indicated 
(filled symbols). PGP (Pretty Good Privacy encrypted communication network) [34:: N = 10680, 
(k) ~ 4.55; P2P (Gnutella peer-to-peer network) [35J: N = 10876, (k) ~ 7.35; PGR (power grid) 
(36], N = 4941, (k) ~ 2.67; APC (astrophysics collaboration) [31J: N = 16706, (k) ~ 14.5. Besides 
the different details of real degree sequences, we can interpret the main features of these nets in 
terms of those observed for the Poisson and power-law classes. 

for similar parameters (k) and N of the Poisson case. In fact its degree distribution decays 
exponentially, (ii) P2P presents power-law decay of the degree distributions, for k > 10, 
with exponent close to 7 = 4. Both plots of P2P are in agreement with those obtained 
for the 7 = 4 class with similar values of N and (k), despite the distributions only share 
in common the power-law tail, (iii) PGP (of size similar to P2P) has a degree distribution 
that decays with exponent 7 < 3 for A; < 50 and 7 ~ 4 beyond [31]. The plot for 1/E vs r 
presents larger values of 1/E than P2P consistent with its (k). However, the plot C vs r of 
PGP deviates from the standard behavior, presenting larger values of C that increase with 
r in a single regime. The interval of allowed values of r is sensitive also to features other 
than the tails. These deviations can be attributed to different initial power law regimes, 
(iv) Finally, APC has a power-law decay with exponent 7 ~ 1 and exponential cut-off for 
k > 50 [31 J . The low and constant plot of 1/E vs r is expected for a network with large 
(k), almost independently of the class of degree distribution. The large values of C are also 
consistent with heterogeneous distributions with large (k). 

Then, despite the different details of real degree distributions, the main observed features 
can be interpreted in terms of the analyzed classes with corresponding values of parameters 
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(k) and N. 



IV. DISCUSSION AND FINAL REMARKS 



For all classes of networks considered, C increases with r in the whole allowed range of r. 
This behavior has already been observed in Ref. |15j . where only positive values of r were 
analyzed and also in Ref. [18] although different definitions of clustering and correlation were 
used. However, we observed that, in the r-ensemble, a non-vanishing clustering coefficient 
C is typically due to finite size effects, such that, in the large size limit, network transitivity 
vanishes as 1/N. In contrast, for power-law networks characterized by r above a threshold, 
a significantly non-null transitivity arises, apparently persisting for large N. 

In any case, since rewiring in the r-ensemble turns C typically small, transitivity does 
not seem to contribute for attaining the prescribed value of r. To identify the factors that 
contribute to r, it is useful to rewrite Eq. (fxl) _ Recalling that (k n ) e = (k n+1 ) / (k) [37], where 
(• • • ) (without subindex) means computed over the degree distribution P(k), then Eq. ([I]) 
becomes 

(k) 2 (kk') e - (k 2 ) 2 
(k)(k 3 ) - (k 2 ) 2 ' { ) 

Following the decomposition made by Estrada [TTJ, notice that the quantity P 3 = 
Y2(k k')(k ~ — 1)' where the sum is performed over all the different pairs of neighboring 
vertices, is the number of paths of length three, then P3 = 3ua + P3, where P3 is the number 
of nontriangular paths of length three (involving four vertices). As done in Eq. ^ for 3ri/\, 
let us scale also P 3 by the number of wedges (paths of length two) P 2 = \ Si h(h — 1), 
defining B = P3/P2 (scaled branching) [T7]. Then Eq. (pi) can be written as 



r 



(k)((k 2 )-(k))(B + l-^ + c) 



(9) 



(k)(k 3 ) - {k 2 } 2 

Expression ^ is determined by the three first raw moments of P(k), and also by B and C 
that are the quantities embodying the information on the linear degree-degree correlations. 

For the Poisson distributed networks, by taking into account the analytical expressions 
for the moments of P(k), it is straightforward to see that 



r = B - (k) + C . 



(10) 
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Clearly, dissassortative correlations are favored by vanishing C. Only for positive r the 
growth of C is convenient, but B — (k) can vary in a wider interval than C (twice wider 
in this case). The existence of two regimes in the increase of C(r), observed in Fig. [TJ is 
consistent with this picture. In other words, Eq. (10) indicates that, in the rewiring process 
of a Poisson network to attain r*, as soon as P(k) is conserved and C remains very small, 
then r is ruled predominantly by B. 

For the power-law distributed networks, some qualitatively similar effects occur, as far as 
the relation of r with B and C is always linear and C is constrained to a narrower interval 
than B. The formation of triangles also in this case contributes only for assortative values 
of r (above the crossover), with values of C larger than in the Poisson case but still small. 
Then also in this case the increase of the branching must drive rewiring. In contrast to 
the Poisson case, the other terms in Eq. ([9]), related to the moments of P(k), might have 
a crucial influence on r because of the divergence, in the infinite network limit, of the nth 
moments for 7 > n + 1. 

Let us analyze the large k max (hence N) limit, setting aside the marginal (logarithmic) 
cases. Here, we use the result k max ~ N. However, if used instead k c ~ jy 1 ^ 7-1 ), the 
conclusions would remain the same. Considering the expressions for the moments (e.g., 
Eq. @), one has, for 3 < 7 < 4: r ~ [B - 0(1)]/ '0(N^), while for 2 < 7 < 3: 
r ~ [B — 0(N 3 ^' y )]/0(N), meaning 0(x a ) ~ ax a , with a > 0. To approach the lower 
limit r = — 1, one must have minimal B. If it is of order greater than the other terms in the 
numerator, then one can not have negative r, because B is non-negative and it will dominate 
the numerator. Thus, negative r can arise only if B is of the same or lower order. But in 
that case r — > in the large N limit. This explains why the lower bound r min tends to 
when 7 < 4 (see Fig. |2](b)). Along this line, however, r min is not expected to vanish when 
7 > 4, but to tend to a small finite value. Similarly, to attain a non-null upper bound of r, 
B needs to grow like the denominator, otherwise, the upper bound will be negative and also 
vanish when iV — > 00, leading to the collapse of the upper bound too. However, this does 
not necessarily happens if B is driven to grow enough during rewiring, which is what seems 
to be happen according to Fig. |2](b). 

The connection between r and distance measures is not so direct analytically. Numerical 
results showed that, for networks with localized distribution of links, changing r modifies 
significantly the mean path length only when correlations are assortative (r > 0.5) and (k) 
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small. These changes could be related to the induced fragmentation, that diminishes by 
increasing (k). Then, the impact of r becomes less important as (k) increases. Meanwhile, 
the influence on the diameter is more dramatic. In power-law networks, the modification 
of the mean path length by r is a bit more marked even if fragmentation is absent for 
(k) > 4, while the diameter is not largely affected. In both cases, the modification of 
characteristic lengths that occur when varying r may affect transport processes and should 
be taken into account either when interpreting or designing numerical experiments on top 
of these networks. 
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